Parametrizations of triaxial deformation and E2 transitions of the wobbling band 
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By the very definition tlie triaxial deformation parameter 7 is related to the expectation values 
of the K = and K = 2 components of the intrinsic quadrupole tensor operator. On the other 
hand, using the same symbol "7", various different parametrizations of triaxial deformation have 
been employed, which are suitable for various types of the mean-field potentials. It is pointed 
out that the values of various "7" are quite different for the same actual triaxial deformation, 
especially for the large deformation; for example, the difference can be almost a factor two for the 
case of the triaxial superdeformed bands recently observed in the Hf and Lu nuclei. In our previous 
work, we have studied the wobbling band in Lu nuclei by using the microscopic framework of the 
cranked Nilsson mean-field and the random phase approximation, where the most serious problem 
is that the calculated B{E2) value is about factor two smaller. It is shown that the origin of this 
underestimation can be mainly attributed to the small triaxial deformation; if is used the same 
triaxial deformation as in the analysis of the particle-rotor model, the calculated B{E2) increases 
and gives correct magnitude compared with the experimental data. 

PACS numbers: 21.10.Re, 21.60.Jz, 23.20.Lv, 27.70.-fq 



I. INTRODUCTION 

The triaxial deformation of atomic nuclei has been one 
of the long standing issues in nuclear structure physics. 
If not going back to the Davydov-Phillipov model [1], the 
triaxial rotor model is very useful to study the spectra 
of odd transitional nuclei [2] . More recently, the model- 
independent sum-rule method [3] is used to analyse the 
Coulomb excitations for the quadrupole collective mo- 
tions, where the triple E2 matrix elements is related the 
the triaxiality parameter 7. Many of Coulomb excita- 
tion measurements revealed clearly non-axial deforma- 
tion [4], although it was difficult to distinguish that the 
non-axiality is of static or of dynamic nature. At high- 
spin states, it was expected that the effect of triaxial 
deformation appears more explicitly in the excitation en- 
ergies and/or M1/E2 transition probabilities, see [5, 6] 
and references therein; for example, the so-called signa- 
ture staggering in odd or odd-odd nuclei was supposed 
to be a good indicator of triaxiality, but the result was 
not so definite as expected. 

However, the situation has been changed quite re- 
cently: The nuclear wobbling motion [7] has been iden- 
tified in some Lu isotopes, ^^^l^ [8-11], ^'^^Lu [12], 
i^'^Lu [13], and ^'^^Lu [14], and it is now believed that 
we are able to study how a triaxially deformed nucleus 
rotates. In fact it has been predicted that the strong 
triaxial deformation appears in this mass region [15-17], 
i.e. the Hf and Lu region, where the wobbling phonon 
excitations have been measured. The triaxial deforma- 
tion predicted in such nuclei is the so-called positive 7 
shape in the Lund convention [18], i.e. nuclei rotate 
about the shortest axis, with the Nilsson deformation 
parameters (£2,7) ~ (0.43,20°), and the associated ro- 
tational sequence is called the triaxial superdeformed 
(TSD) band [19]. 



In order to pin down how much the TSD nucleus de- 
forms triaxially, in which the wobbling excitation is mea- 
sured, one needs to have recourse to some model. The 
standard one is the triaxial rotor model [1], where the 
high-spin states should be considered [7] . In this model, 
the effect of triaxial deformation appears primarily in the 
three different moments of inertia and the two intrinsic 
quadrupole moments of nucleus. The formers are mainly 
responsible to the excitation energies, while the latters 
are directly reflected in the E2 transition probabilities. 
Compared with the moments of inertia which are sensi- 
tive to other factors, like the pairing correlations, the E2 
transition probabilities are more direct and robust quan- 
tities, and therefore, the measurements oi B{E2) [20, 21], 
not only that of in-band (intraband) but of out-of-band 
(interband) between the yrast TSD band and the ex- 
cited wobbling band, are crucial to obtain the informa- 
tion about the triaxial deformation. The detailed study 
by using the particle-rotor model [8, 22, 23], where an 
odd «i3/2 proton is coupled to the triaxial rotor, which 
is suitable for the description of odd Lu TSD bands, re- 
vealed that the observed ratio of out-of-band to in-band 
B{E2ys, B{E2)out/B{E2)in, is consistent to the triaxi- 
ality parameter 7 w 20°. 

In this respect, it is worthwhile mentioning Ref. [24], 
where it is discussed that the out-of-band B{E2) from the 
one-phonon wobbling band is not enough to distinguish 
the 7-soft and 7-rigid triaxialities in nature. Further, it 
is proposed that the measurement of the B{E2) between 
the odd-spin and even-spin members of the wobbling ex- 
citations is crucial to distinguish the two. Here even- 
even nuclei are considered in [24] so that the observed 
one-phonon wobbling band corresponds to the odd-spin 
members for the even-spin yrast TSD band. In some Lu 
isotopes the so-called two-phonon wobbling bands have 
been observed [11, 12], which corresponds to the even- 
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spin members of the wobbling excitations for even-even 
nuclei. The B{E2) values from the two-phonon to one- 
phonon wobbling band are measured to be about two 
times the values from the one-phonon wobbling to yrast 
band, which clearly fits the picture of the 7-rigid model 
rather than the 7-soft model. Although the negative 7 
shape is assumed in Ref. [24], which is believed to be op- 
posite to what is measured in Lu isotopes (see below), 
it does not affect the main issue; whether the 7-soft or 
7-rigid model is superior. In this way, it is more likely 
that the triaxial deformation is of static nature, i.e. the 
rigid model is more suitable, and it is meaningful to ask 
how much triaxial deformation the observed TSD bands 
have. 

It should, however, be noticed that the conventional 
macroscopic rotor model with irrotational moments of in- 
ertia has an essential problem; the rotor rotates around 
the intermediate axis, which corresponds to the nega- 
tive 7 shape, and conflicts to the measured B{E2) ratio. 
Therefore, in Refs. [8, 22, 23], the largest and intermedi- 
ate moments of inertia are interchanged to simulate the 
positive 7 shape. It is well-known that the nuclear mo- 
ments of inertia are neither that of the rigid-body nor of 
the irrotational fluid, and microscopic models are neces- 
sary for their proper description. Thus, we have investi- 
gated the wobbling motion in the Lu region [25 27] by 
employing a microscopic framework, the cranked mean- 
field and the random phase approximation (RPA) [28- 
30]. This approach is suitable to describe the vibra- 
tional excitations in the rapidly rotating nuclei, see e.g. 
Refs. [31-35]; and have been used to study possible wob- 
bling excitations in normal deformed nuclei in our previ- 
ous works [36, 37], and more recently in Refs. [38, 39]. 

In our approach the residual interaction is chosen con- 
sistently to a given mean-field, and there is no ambiguity 
of the force parameters for the RPA calculation. It has 
been found [25-27] that the RPA solutions, which can 
be nicely interpreted as the wobbling phonons, do ex- 
ist in the Lu region for a suitable range of deformation 
parameters corresponding to the prediction of the TSD 
bands. The calculated excitation energies is in a rea- 
sonable range in comparison with the experimental data, 
though they do not fit the data very precisely. However, 
the calculated B{E2) ratios are systematically smaller by 
about a factor two to three, as long as is used the triaxial- 
ity parameter 7 w 20°, which is predicted by the Nilsson- 
Strutinsky calculations [17, 19] and is used in the anal- 
ysis of the B(E2) ratio by the particle- rotor model [23]. 
The main issue of the present work is to discuss why 
our RPA calculations of the wobbling motion underpre- 
dict the B{E2) ratio. In the course of discussion, it is 
clarified that the various definitions of the triaxiality pa- 
rameter 7 give rather different values and one has to be 
very careful when talking about the triaxial deformation 
especially for a larger deformation like in the case of TSD 
bands. 

The paper is organized as follows: Various existing def- 
initions of the triaxiality parameter 7 are reviewed, and 



their values for a given shape are calculated and com- 
pared in §11. Some general discussion on the relations 
between various definitions are also given there. After 
discussing the difference between the 7 values used in 
the Nilsson-Strutinsky calculations and the particle-rotor 
model, in §111, it is shown that our RPA calculation gives 
a correct magnitude of the B{E2) ratios if is used the 
corresponding triaxial deformation to the analysis of the 
rotor model. §IV is devoted to the summary. A part of 
the present work was presented in some conference re- 
ports [40, 41]. 

II. PARAMETRIZATIONS OF TRIAXIAL 
DEFORMATION 

The amount of triaxial deformation is usually desig- 
nated by the triaxiality parameter 7, but there are var- 
ious definitions for it. In this section, we discuss the 
relations between them and show how large their dif- 
ferences are for a given shape. It should be mentioned 
that this problem has been already discussed in Ref. [42] 
(Appendix B) for the volume- conserving anisotropic har- 
monic oscillator potential. The present study generalizes 
its conclusion to more realistic potentials. 

A. Basic definition based on the quadrupole 
moments 

One of the most important characteristics of nuclei 
with (static) triaxial deformation is the existence of two 
distinct intrinsic quadrupole moments. In this paper we 

assume the intrinsic z-axis as a qiiantization axis and the 
X-axis as a rotation axis, and define the two moments [7], 

^ ^ r r W 
_ Q2 ^ ^'^{Q22) = ^Jlj {x' - y')p{r)d'r, 

where Q2K {K = 0, ±1, ±2) are the usual quadrupole op- 
erator in the intrinsic frame of the deformed nucleus, and 
p{r) is the nucleonic density. These two moments are di- 
rectly related to the in-band and out-of-band B{E2) val- 
ues of the wobbling band according to the rotor model [7] , 
and their measurements in the Lu isotopes [20, 21] 
uniquely determine these moments. In place of the two 
moments, equivalent two quantities, the magnitude of 
moments Q and the triaxiality parameter 7 are usually 
used: 

(5o = Qcos7, -V2Q2 = Qsm.'y. (2) 

Here we follow the Lund convention [18] of the sign of 
the triaxiality parameter 7, which is opposite to that 
of Ref. [7]. In the following we are mainly concerned 
with the parameter 7, and consider the in-band and out- 
of-band B{E2) ratio of the one-phonon wobbling band. 
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which is independent of the magnitude of the moment Q; 
see the next section. 

The triaxiaUty parameter defined above reflects the nu- 
clear density distribution and we call it "7(den)" in this 
work; i.e. 



tan 7 (den) 



V2(Q22) 

(Q20) 



(3) 



Since 7 (den) depends on the calculated single-particle 
wave functions of the average nuclear potential, and on 
the configuration of each nucleus, it is more convenient 
to introduce an another parameter "7(geo)", which is 
more directly related to the geometric shape of nucleus 
mathematically defined by the two dimensional surface 
E: 



tan7(geo) = — 



(Q20)un 



(4) 



where ( )u„i means that the expectation value is taken 
with respect to the sharp-cut uniform density distribu- 
tion. 



Puni(r) 



Po for r inside the surface S, 
otherwise. 



(5) 



coincide typically within 10% except in the small 7 re- 
gion, where both 7 (den) and 7(geo) become zero and 
the ratio is numerically unstable. This agreement cor- 
responds to the so-called shape consistency between the 
density and the potential [43] , which has been tested both 
for the Nilsson and Woods-Saxon potentials for the axi- 
ally symmetric deformations [44, 45]. Therefore, we can 
practically use 7(geo) in place of 7(den). 

The definition of triaxiality parameter, 7(den) or 
7(geo), is basic or fundamental in the sense that it is 
directly related to the B{E2) values of the triaxial rotor. 
In the microscopic calculations by means of the Hartrcc- 
Fock or Hartree-Fock-Bogoliubov method, 7 = 7(den) is 
the only definition of the triaxial deformation, and no 
confusion exists. It is, however, quite often that one 
starts from some average nuclear potential, which is tri- 
axially deformed, and calculates the potential energy sur- 
face or the total routhian surface to determine the self- 
consistent deformation as in the case of the Strutinsky 
method. In such cases, a different type of definitions of 
the triaxiality parameter 7 has been used, which is cho- 
sen conveniently to specify the triaxial deformation of the 
potential. We call this third type of definition "7(pot)": 
We consider two well-known definitions in the foUowings, 
depending on the employed potential. 
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FIG. 1: The ratios 7(dcn)/7(geo) arc plotted as functions of 
7(geo) for the cases of the Nilsson (solid line) and Woods- 
Saxon potentials (dashed line). The calculations have been 
performed for the TSD band in ^^^Lu, and the parameters 
used for the Nilsson potential are €2 = 0.43, e4 = 0.0, and 
^n,p = 0.3 MeV, and those for the Woods-Saxon potential 
are /32 = 0.42, = 0.034, and Z\n,p = 0.3 MeV, which has 
almost the same shape as the case of the Nilsson potential 
at 7(geo) « 10° in the minimum of the Nilsson-Strutinsky 
calculation. 

If the density is calculated within the mean-field ap- 
proximation by using the average single-particle poten- 
tial which has the same shape as the uniform density (5), 
then these two triaxiality parameters, 7(den) and 7(geo), 
agree very well. In Fig. 1 the ratio 7(den)/7(geo) is 
shown as a function of 7(geo), and we can see that they 



B. 7(pot) defined in the Nilsson potential 

As a definite example of 7(pot), we take the Nils- 
son potential, or the modified oscillator potential, as 
a mean-field potential, i.e. 7(pot:Nils). The deforma- 
tion parameters in the Nilsson potential [18, 46, 47] 
considered in the present work are (£2,7,64), which de- 
fine the deformation of the velocity independent part of 
po tential through the si ngle-str etched coordinate, r' = 
("v/wcn/woa;, ^July/uJoy, ^Juiz/uJoz), as 

V{r) = iMa;oa;v(e2,7,e4)r'^x 

(1- ^ C2kY2k{^')- C4ifF4i^(J^')), 

K=0.±2 K=0.±2,±4, 

(6) 

where ojq is the frequency of the spherical potential, 
'^v(e2, 7, £4) is determined by the volume conserving con- 
dition, ri' is the solid-angle of coordinate r', and the 
coefficients c's are given by 

C20 = £2 cos 7, 

C22 = C2-2 = £2 sin 7, 

' C4o = ^e4(5cos2 7 + 1), 

C42 = C4_2 = - ^^^O"" 64 cos 7 sin 7, 
^ C44 = C4-4 = £4 sin^ 7. 

Note that the three frequencies, cjj,, ujy, and cj^, are 
given by Eq. (9) below. The nuclear shape S in this 
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case is defined as an equi-potcntial surface of the po- 
tential, V{r) — const., and is uniquely determined once 
the parameters (e2, 7, 64) with 7 — 7(pot:Nils) are given. 
It is straightforward but rather complicated because of 
the use of the single-stretched coordinate r' in prac- 
tice. Using this potential either the triaxiality parameter 
7 (den) or 7(geo) defined in the previous subsection can 
be calculated as functions of these potential parameters 
(e2,7 = 7(pot:Nils),e4). 




10 20 30 40 50 60 

7(pot:Nils) 

FIG. 2: The triaxiality parameters 7(den) (dashed line) and 
7(geo) (solid line) are shown as functions of the 7 parameter 
in the Nilsson potential, 7(pot:Nils) (the thin diagonal line 
is just for a guide to the eyes). The calculations have been 
performed for the TSD band in ^^"^Lu, and the parameters 
used are ei = 0.43, £4 = 0.0, and Z\n,p = 0.3 MeV. 




10 20 30 40 50 60 

7(pot:Nils) 

FIG. 3: The triaxiality parameters 7(geo) with fixed £2 and £4 
deformation parameters are shown as functions of 7(pot:Nils) 
(the thin diagonal line is just for a guide to the eyes). The 
dotted line is for the case with £2 = 0.2, the dashed with 
£2 = 0.4, and the solid with £2 = 0.6, respectively. The £4 
parameter is set for all the cases. 

In Fig. 2, are depicted the relation between 7 (den) 

and 7(pot:Nils) and that between 7(geo) and 7(pot:Nils) 
at a given (€2,64), which is suitable for the TSD band 



in °^Lu. As is already shown in Fig. 1, 7(den) and 
7(geo) are very similar, but 7(pot:Nils) is quite different: 
7(pot:Nils) = 20° corresponds to 7(den) Ri 7(geo) ~ 11°, 
so that the difference can be as much as about a factor 
two. In order to see how the difference between 7(geo) 
and 7 (pot: Nils) changes for different £2 values, three cases 
with €2 = 0.2, 0.4, 0.6 are shown in Fig. 3. It is clear that 
the difference becomes larger for larger £2 deformations: 
7(geo) is only about 10° even though 7(pot:Nils) is put 
30° in the case of the superdeformed band, £2 « 0.6. 

In the case of £4 = 0, the Nilsson potential reduces to 
the anisotropic harmonic oscillator potential except for 
the and I ■ s terms, which are irrelevant for the defi- 
nition of the nuclear shape. It is instructive to consider 
such a case in order to understand the difference shown 
in Fig. 3. Then the shape is a volume-conserving ellipsoid 
defined by simple equations. 
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2 2 



const., with JJc 



(8) 



The frequencies uji [i = 1,2,3) for the x, y, z-directions, 
which are inversely proportional to the lengths of the 
ellipsoid along these axes, are given by 



(7 + Y*)^ ' ^ = 7(pot:Nils). 



iOi = (jO^ \ ~ 3^^ 

(9) 

Therefore 7(geo) and 7(pot:Nils) are related through 
(1 = 1,2,3); 



tan7(geo) = ^ _2 



tan7(pot:Nils) = 



-2 ' 



(10) 



(11) 



which relate these two 7's for a given value of £2. In the 
limit of small deformation parameters, £2,|7| <C 1, it is 
easy to confirm 



7(geo) ~ ( 1 - ^£2 ) 7(pot:Nils). 



(12) 



Namely, the slope of curves at the origin in Fig. 3 changes 
with £2 with a rather large factor |, and this clearly ex- 
plains that 7(geo) is smaller than 7(pot:Nils) more than 
a factor two when £2 is as large as 0.4 like in the case of 
the TSD band. 

In order to see how these different definitions of two tri- 
axiality parameters, 7(geo) and 7(pot:Nils), change the 
appearance of potential energy surface, we show an exam- 
ple in Fig. 4. Here the £4 parameter is chosen to minimize 
the potential energy at each (£2,7) mesh points. The pa- 
rameter 7(gco) depends not only on (£2, 7(pot:Nils)) but 
on £4, and it is impossible to calculate the (£2,7(geo)) 
mesh points before the minimization with respect to £4. 
Therefore, we made an approximation to set £4 = 
when we prepare the (£2,7(geo)) mesh points from the 
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FIG. 4: Potential energy surface obtained by the cranked Nils- 
son Strutinsky calculation for the {n,a) = -1-1/2) config- 
uration in ^®^Lu at 7 = 41/2+. The energy between contours 
is 250 keV. The triaxiality parameter 7 = 7(pot:Nils) is used 
as usual in the upper panel, while the triaxiality parameter 
7 = 7(geo) in the lower panel. 



(£2, 7(pot:Nils)) mesh points. As is clear from the fig- 
ure, the surface is squeezed to tlic 7 = axis at larger 
deformation, and apparently the TSD minimum moves 
to smaller triaxial values. In Fig. 4, only the 7 param- 
eter is replaced from 7(pot:Nils) to 7(geo). However, it 
may be better to replace £2 to the other parameter corre- 
sponding to the magnitude Q in Eq.(2) in order to make 
the meaning of the quadrupole deformation clearer. Con- 
straint Hartree-Fock(-Bogoliubov) type calculations are 
necessary for such a purpose. Then the method becomes 
much more involved and the simplicity of the Strutinsky 
type calculation may be lost. 

It may be worthwhile mentioning, here, the specific 



model composed of the spherical Nilsson potential and 
the QQ force as an effective interaction. The velocity- 
independent part of the Hartree potential in this model 
is given by 



V{r) = -Mujq r - a2oQ20 - ol22{Q22 + Q2-2), 



(13) 



where the Hartree condition requires a2K = X {Q2k) 
{K = 0,2) with X being the QQ force strength, and 
the potential reduces also to the anisotropic harmonic 
oscillator (but the volume conservation condition is not 
necessarily satisfied in this Hartree procedure). If the 
model space is not restricted, a2K is proportional to the 
two intrinsic quadrupole moments in Eq. (1), and then, 
the triaxiality parameter of the density type (3) in this 
case is related to — a/2 a22/Q;2o, which can be expressed 
in terms of (i = 1, 2, 3) as 



tan7(den:QQ) 



2u;2_^2 



,,2 ■ 



(14) 



It is usual to parametrize the potential by the deforma- 
tion parameter (5, 7) with 7 = 7(den:QQ)), in terms of 
which the frequencies are written in the form, 



^t]], 7 = 7(den:QQ). 



1 — —5 cos 

(15) 

Thus, the triaxiality parameters 7(geo) and 7(den:QQ) 
are quite different in this case: In the limit of small de- 
formation parameters, ^,|7| <C 1, they are related like 



7(geo)«(l-2<5)7(den:QQ). 



(16) 



Namely, the difference is even larger than that between 
7(geo) and 7(pot:Nils). This means that the shape con- 
sistency between the density and the potential is strongly 
violated in the QQ force model, if the full model space is 
used in the Hartree procedure. Considering this on top 
of the fact that the volume conservation condition is not 
guaranteed, the model is not realistic at all. 

Actually, the QQ force model is supposed to be a model 
in a restricted model space, and the contributions of the 
"core" should be added to {Q2k) calculated within the 
model space in order to obtain the intrinsic moments 
of the whole system. Then the problem of breaking 
the shape consistency may not be the real problem. In 
fact, the core contributions of the quadrupole operators, 
i.e. the expectation values for the closed shell configu- 
rations in the anisotropic harmonic oscillator potential, 
have different dependence on the deformation; it can be 
shown [42] that they lead to the triaxiality parameter 
7 = 7(den:core), which satisfies 



tan 7 (den: core) 



\/3(u 



-1 ■ 



(17) 



In Ref. [42] (see Appendix B of this reference for de- 
tails, but note the different notations used there), these 
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differences between various types of the deformation pa- 
rameters defined in tlie liarmonic osciilator potential 
were already discussed, where not only the parameter 
7 but also the other one of a pair of the parameters, 

^ ~ '\/^^/(Sfc=i ''^)' wsre considered. It was already 
pointed out that the triaxiality parameters for a given 
shape in various types of definition take quite different 
values. 



C. 7(pot) defined in tiie Woods-Saxon potential 

As an another example of 7(pot), the parametrization 
of deformation in the Woods-Saxon potential is consid- 
ered, i.e. 7(pot:WS). Actually, it is not restricted to 
the Woods-Saxon potential, but is more general as one 
can see in the following. The deformed Woods-Saxon 
potential considered in this work is parametrized by the 
deformation parameters, (/32,7, /34), and defined [48-50] 
by 



V{r) 



1 + exp(dists(T')/a) 



(18) 



where distE(T') is the distance between a given point r 
and the nuclear surface S, with a minus sign if r is inside 
S, which is defined by the usual radius to solid-angle 

relation, r = R{il); 

R{n) = R,{^2n,Mx 

(l+ ^ a2KY2K{n)+ aiKY4K{^)), 



K=0,±2.±4 



(19) 

where i?v(/?2, 7, Pa) is determined by the volume conserv- 
ing condition, and the coefficients a's are given by 



020 = P2 COS 7, 

a22 = a2-2 = - /32 sin 7, 
^40 = ^/?4(5cos^7-|- 1), 



■1' 

y|/34C0S7sin7, 
a44 = a4_4 = ^J^Pa sin^ 7- 



a42 



6 ' 
04-2 



(20) 



Apparently the surface E is given as an equi-potential 

surface at the half depth, V{r) = ■^Vq, and it is directly 
related to (/32,7,/?4) with 7 — 7(potWS). 

As in the case of the Nilsson potential, the triaxial- 
ity parameter 7 (den) and 7(geo) are shown as functions 
of 7(pot:WS) in Fig. 5, which are suitable to the TSD 
band in -"^^^Lu, just in the same way as in the case of the 
Nilsson potential in Fig. 2. Again, 7 (den) and 7(geo) are 
very similar, but they are quite different from 7(pot:WS): 
7(den) w 7(geo) « 13° when7(potWS) = 20°. In Fig. 6, 
the relation between 7(geo) and 7(pot:WS) at three dif- 
ferent cases of {/32,04) deformations, corresponding to 
Fig. 3, are also depicted. Although the differences be- 
tween 7(geo) and 7(pot:WS) are not so dramatic as those 
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FIG. 5: The triaxiality parameters 7(den) (dashed line) and 
7(geo) (solid line) are shown as functions of the 7 parameter 
in the Woods-Saxon potential, 7(pot:WS) (the thin diagonal 
line is just for a guide to the eyes). The calculations have been 
performed for the TSD band in ^^''Lu, and the parameters 
used are P2 = 0.42, 04 = 0.034, and An,p = 0.3 MeV; these 
deformation parameters gives almost the same shape as that 
used in the Nilsson potential in Fig. 2 at 7(geo) « 10°. 
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FIG. 6: The triaxiality parameters 7(geo) with fixed 13-2 and /34 
deformation parameters are shown as functions of 7(potWS) 
(the thin diagonal line is just for a guide to the eyes). The 
dotted Une is for the case with /32 = 0.217 and /34 = 0.017, 
the dashed with (32 = 0.445 and (Sa = 0.075, and the soUd 
with (32 = 0.685 and B4, = 0.190, respectively. These sots 
of parameters almost correspond to the cases of t2 = 0.2, 
£2 = 0.4, and £2 = 0.6 with £4=0 used for the Nilsson 
potential in Fig. 3 at 7(geo) = 0. 



between 7(geo) and 7(pot:Nils), they are still consider- 
ably large. In the case of the parametrization of the 
nuclear surface in Eq. (19), ((32if)uni can be easily calcu- 
lated; 

((52i^)uni = \j R{nfY2K{^)dn. (21) 

Then, it is straightforward to see, in the small deforma- 



7 



tion limit, /32,|7| ^ 1, with /34 = 0, that 

7(geo) « (l - 7(pot:WS). (22) 



Taking into account the relation, /32 ~ y £2 in the 
small deformation limit, the proportionality constant in 

front of £2 corresponds, in this case, to \/ ^P- x < / = 



497r Y 45 

I « 1.14, which is smaller than | = 1.5 in Eq. (12) 
for the Nilsson potential, but is still appreciably large. 
This explains qualitatively the increase of the difference 
between 7(geo) and 7(pot:WS) for larger deformations 
as is shown in Fig. 6. 



III. B{E2) RATIO OF THE WOBBLING BAND 

As it is discussed in the previous section, the two in- 
trinsic quadrupole moments should be determined in or- 
der to deduce the triaxial deformation. In the case of 
the wobbling excitations, it is enough to measure the two 
B{E2ys, B{E2)in and S(S2)out; the AI = -2 intraband 
E2 transitions within the wobbling band and AI = ±1 
interband E2 transitions from the one-phonon wobbling 
band to the yrast TSD band, respectively. According 
to the rotor model [7], the magnitude of the moment Q 
in Eq.(2) is factored out in the two B{E2ys and their 
ratio, B{E2)out/B{E2)in, is directly related to the triax- 
iality parameter 7(den). This B{E2) ratio is straightfor- 
ward to measure from the experimental point of view; it 
is given directly by the 7-ray branching ratio if the in- 
formation of the mixing ratio is provided. In contrast, 
the life time measurement is necessary to obtain B{E2) 
values themselves, which is not an easy task generally. 
Although the life time measurements have been done in 
some TSD bands [20, 21] recently, so that we can study 
both B{E2)out and B{E2)in separately, we concentrate 
upon the B{E2) ratio in the present work. 

In Fig. 7, the experimental B{E2) ratio of the one- 
phonon wobbling band in ^^^Lu [21] is compared with the 
results of the particle-rotor model calculation in Ref . [23] . 
Important parameters of the model are three moments of 
inertia, J^x^ Jy, Jzi and the triaxiality 7 = 7 (den); an 
overall factor of the formers is irrelevant to the B{E2) ra- 
tio and they are fixed to be Jx '■ Jy '■ Jz = 145 : 135 : 50 
taken from [23], while we take five values 7(den) = 10°, 
15°, 20°, 25°, and 30° in order to show the dependence of 
the B(E2) ratio on the 7 values. Other parameters, the 
chemical potential of the odd 113/2 proton, A/k = —1.532, 
and the pairing gap, A/k = 0.3, with k = 3 MeV are 
also taken from [23]. The calculated ratios are mono- 
tonically decreasing functions of spin if all the model pa- 
rameters arc held fixed. This decrease is characteristic in 
the rotor model, see Eq. (23) below. Although the spin- 
dependence is somewhat different, the average value of 
the B{E2) ratio can be reproduced if we take the value 
7 « 20°. Therefore the triaxial deformation of ^®^Lu is 
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FIG. 7: Tiic B{E2) ratio, B{E2:I I - i)out/B{E2: 1 
I — 2)in, of the wobbling band in ^^'^Lu. The experimental 
data [21] is compared with the calculations by the particle- 
rotor model [23]. Five dotted curves are the results with the 
triaxiality parameter 7(den) = 10°, 15°, 20°, 25°, and 30° 
from the lowest to the highest ones, respectively. 



deduced to be 7 « 20°, which is one of the main con- 
clusions of Ref. [23] . This conclusion remains valid even 
if the parameters of the moments of inertia are changed 
in a reasonable range. While the existence of the odd 
proton brings about important corrections to the energy 
spectra, its effect on the B{E2) is very small [23, 51]. 
Thus, the wobbling phonon treatment of the simple ro- 
tor model in Ref. [7] gives a good approximation to the 
B{E2) ratio at high-spin states, leading to the following 
expression; 

B(£;2:/^/±l)out 



B{E2: 1 ^ I - 2)in 

2 f Wz sin(7 -I- 60°) ^ Wy sin 7 
/V ^WyWz cos(7 -I- 30°) 



7 = 7(den), 



(23) 

where the quantities Wy , are related to the three mo- 
ments of inertia through 



^{jx/jy-iy^. 



(24) 



Note that the 7 — > 7 -|- 1 transitions are quenched for 

the positive 7 shape, and in fact only the / ^ / — 1 
transitions are observed in experiment in the Lu nuclei. 
It is easy to check that the results of calculations in Fig. 7 
can be understood nicely by this simple expression; the 
decrease of the B{E2) ratio as a function of spin is due 
to the 1/7 dependence in Eq.(23), and the ratio increases 
quickly as a function of the triaxiality 7 = 7(den). 

It may be interesting to note that the measured B(E2) 
ratio is almost constant as a function of spin or even 
increases at highest spins, which is quite different from 
that of the rotor model calculation. It indicates that 
the parameters of the model are changing as the spin 7 
increases. Since the ratio is sensitive to the triaxiality pa- 
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ramctcr 7 = 7 (den), it is natural to consider that 7 (den) 
is spin dependent and increases with spins [23] among 
others. A preliminary investigation for such a possibility 
has been reported in Ref. [41], where is employed a micro- 
scopic framework, the cranked Woods-Saxon mean-field 
and the random phase approximation (RPA). 

Next let us turn to the discussion on our microscopic 
calculations in Refs. [25 27], which are based on the 
cranked Nilsson mean-field and the random phase ap- 
proximation (RPA). In these calculations the triaxial- 
ity parameter 7 = 20° was employed but the resultant 
B{E2) ratios were too small by a factor two to three. We 
were wondering about possible reasons; does the results 
of RPA calculation deviate from those of the rotor model 
so much? However, it has been shown in Ref. [52] that 
the RPA calculation reproduces the result of the rotor 
model rather well in the case of the precession bands, 
which are nothing but the rotational bands built upon 
the high-if isomers and can be interpreted as a simi- 
lar motion to the wobbling excitation, where the angu- 
lar momentum vector fluctuates about the main rotation 
axis [53]. Now the reason of the small calculated B{E2) 
ratio is clear from the argument of the previous section: 
The triaxiality 7 used in our calculations is that of the 
Nilsson potential, 7(pot:Nils), in §11 B, while the 7 used 
in the rotor model is 7(den) in §11 A. As is discussed in 
the previous section, the difference between them for the 
same shape is very large for large deformations like in the 
case of the TSD bands; 7(pot:Nils) k, 20° corresponds to 
7(den) « 11°, and in order to perform the same calcu- 
lation as the rotor model with 7 (den) w 20° one has to 
employ 7(pot:Nils) « 30° according to Fig. 2. It should 
also be mentioned that we have used five major oscillator 
shells, A^osc = 3 — 5 for proton and Nose = 4 — 8 for neu- 
tron in the calculation in Refs. [25-27], which were not 
enough; since the ^13/2 proton orbits are occupied in the 
TSD band, the inclusion of A^osc = 8 proton quasiparticle 
states are necessary in the RPA calculational step. 

In Fig. 8 we depict the new results of calculation em- 
ploying 7(pot:Nils) = 20° and 30° with using the full 
model space; all orbits in the oscillator shell A^osc = — 9 
for both protons and neutrons are included. The proce- 
dure and the other parameters in the calculation are the 
same as in the previous work [25]; 62 = 0.43, £4 = 0, 



and the pairing gaps Z\n,] 



0.3 MeV. The result of 



the previous calculation, the dashed line, and that of 
the particle- rotor model with 7 (den) = 20°, the dot- 
ted line, are also included. Our previous calculation is 
smaller than the experimental data partly because of 
the small model space, but its effect is about 20%; the 
large difference is mainly due to the fact that we have 
used 7 = 7(pot:Nils) = 20° in the previous calcula- 
tion, which corresponds to much smaller triaxiality than 
7 (den) — 20° in the particle-rotor model. The result 
with 7(pot:Nils) = 30° almost coincides with that of the 
particle-rotor calculation using 7 (den) = 20°, because 
the values Jx '■ Jy '■ Jz of the microscopically calculated 
moments of inertia accidentally take similar values in the 
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FIG. 8: The B{E'i) ratio, B{E2:I I - i)out/B{E2: 1 
I — 2)in, of the onc-phonon wobbling band in ^^-^Lu. The 
experimental data [21] is compared with the calculations by 
our microscopic RPA approach. The lower solid line is the 
result with 7 = 7(pot:Nils) = 20°, while the upper solid line 
is with 7(pot:Nils) — 30°; the full model space is used for both 
of them. The dashed lino is the result with 7(pot:Nils) = 20° 
and using only the five major oscillator shells, corresponding 
to the previous calculation in Ref. [25]. The dotted line is 
the same as that in Fig. 7, the particle-rotor calculation with 
7 = 7(den) = 20°, depicted for a reference. 
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FIG. 9: The 7 dependence of the microscopically calculated 
B{E2) ratio and the excitation energy of the one-phonon wob- 
bling band in '^^^Lu at spin 7 = 51/2. The solid lines are the 

B(E2) ratios and the dotted lines arc the energy in McV. The 
horizontal solid and dotted lines designate the experimental 
values. 



relevant spin range [25]. In order to see the 7 dependence 
of the results, we show, in Fig. 9, the B{E2) ratio and 
the excitation energy of the one-phonon wobbling band 
in ^^^Lu at spin / = 51/2 as fmictions of the triax;ial- 
ity 7 = 7(pot:Nils). Although the excitation energy is 
rather flat in the range, 20° < 7(pot:Nils) < 30°, the 
B{E2) follows the behaviors of Eq. (23) if the relation 
between 7(pot:Nils) and 7(den) is taken into account. 
The excitation energy can be expressed in terms of the 
three moments of inertia in an usual way [7, 30], but 
their 7 dependences are not simple like the irrotational 
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one [37] and lead to a rather weak 7 dependence of the 
excitation energy in this case. Here we have only shown 
the example of the wobbling excitation in ^®^Lu, but we 
have confirmed that these properties of the wobbling-like 
RPA solution are general, and can be applied for other 
cases in the Lu region. 

Thus, our microscopic RPA calculations give more or 
less the same results as that of the macroscopic particle- 
rotor model, if the corresponding magnitude of the tri- 
axial deformation is employed. It shows that the RPA 
calculation of the wobbling excitation in the Lu region 
leads to the behaviors of B{E2)^s that are given by the 
macroscopic rotor model; namely the out-of-band B{E2) 
can be related to the static triaxial deformation. This is 
non trivial since the out-of-band B{E2) is calculated by 
the RPA transition amplitudes of the non-diagonal part 
of the quadrupole operators, Q2i^ and (522^ PS, 30]. It 
has been shown [37] that the RPA wobbling theory of 
Marshalek [30] gives the same expression of the out-of- 
band B{E2) as that of the rotor model, if the RPA wob- 
bling mode is collective enough that the quantity 
defined in Eq. (4.29) in Ref. [37] satisfies c„=wob = 1- 
In the previous calculations [25-27, 37], the employed 
model space was too small [54] to give c„=wob = 1) but 
we have confirmed that c„=wob ~ 1 is satisfied within 1% 
in the present full model space calculations. Recently, 
this criterion, c„=wob ~ 1) has been used to identify the 
wobbling-like solution out of many RPA eigenmodes, and 
shown to be very useful [38]. 

IV. SUMMARY 

In this work, we have first discussed the differences of 
the various definitions of the triaxiality parameter 7. The 
most basic among them is defined through the two intrin- 
sic quadrupole moments, 7 (den) in Eq. (3) for each con- 
figuration of a particular nucleus, or 7(geo) in Eq. (4) for 
a given nuclear shape. It has been found that these two 
coincide in a good approximation (the nuclear shape con- 
sistency). In the Hartree-Fock(-Bogoliubov) type calcu- 
lations, where the nuclear mean-field is determined self- 
consistently by a suitably chosen effective interaction, the 
parameter 7 (den) is the only possible definition of triax- 
ial deformation. However, there is an another type of 
mean-field calculations, i.e. the Strutinsky macroscopic- 
microscopic method, where one starts from a suitably 
chosen average potential, in which the nuclear defor- 
mation is parametrized in various different ways. We 
have considered the two widely adopted ones, the Nilsson 
type and the Woods-Saxon type parametrizations. The 
triaxiality parameters associated with these potentials. 



7(pot:Nils) in Eqs. (6)-(7) and 7(pot:WS) in Eqs. (18)- 
(19), arc compared with the density type 7(den) and 
7(geo). Conspicuous differences between the potential 
type and the density type 7's, e.g. 7(pot:Nils) vs. 7(geo), 
have been found especially for larger deformations, e.g. 
the triaxial superdeformed states. It is also investigated 
how the differences between various definitions come out 
by evaluating their relations explicitly in the small defor- 
mation limit. Therefore we have to be very careful about 
which definition is used in quantitative discussions of the 
triax;ial deformation. 

Next, we have investigated the out-of-band to in-band 

B{E2) ratio of the one-phonon wobbling band, which is 
measured systematically in the Lu region and is sensitive 
to the triaxial deformation. The macroscopic particle- 
rotor model [23] is used to deduce the triaxial defor- 
mation from the experimental B{E2) ratio, which leads 
to 7 = 7(den) « 20° on average. On the other hand, 
we performed the microscopic RPA calculation [25, 26] 
with using 7 = 7(pot:Nils) « 20° corresponding to the 
TSD minima in the cranked Nilsson-Strutinsky calcula- 
tion [17, 19], but we obtained too small B(E2) ratio com- 
pared with the experimental data. It has been found that 
the reason of the underestimation of our previous micro- 
scopic calculation is mainly due to the different triaxial 
deformation used: We have used 7(pot:Nils) « 20° corre- 
sponding roughly to 7(den) w 11°, which is much smaller 
than 7(den) w 20° in the rotor model calculations. If the 
proper triaxiality corresponding to 7(den) « 20° is used, 
our RPA calculation can nicely reproduce the magnitude 
of the measured B{E2) ratio in the same way as in the 
macroscopic particle-rotor model. 

It should, however, be emphasized that an important 
problem remains: The predicted triaxial deformations, 
7 (den) K. 11°, by the cranked Nilsson-Strutinsky calcula- 
tions [17, 19] for the TSD bands in the Hf, Lu region are 
too small to account for the measured B{E2) ratio of the 
wobbling excitations. We believe that this is a challenge 
to the existing microscopic theory. An another thing we 
would like to mention is that the measured B{E2), both 
the out-of-band and in-band B{E2ys, seems to indicate 
that the triaxial deformation is changing as a function of 
spin; it increases at higher spins [23]. We have recently 
developed a new RPA approach [41] based on the Woods- 
Saxon potential as a mean-field, which is believed to be 
more reliable than our previous calculations employing 
the Nilsson potential. The result of calculations and dis- 
cussions including the issue of the change of the triaxial 
deformation suggested by the out-of-band as well as in- 
band B{E2)^s will be reported in a subsequent paper; see 
Ref. [41] for a preliminary report. 
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